Data-driven derivation of the turbulent energy cascade generator 
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Within the framework of random multiplicative energy cascade models of fully developed turbu- 
lence, expressions for two-point correlators and cumulants are derived, taking into account a proper 
conversion from an ultrametric to an Euclidean two-point distance. The comparison with two-point 
statistics of the surrogate energy dissipation, extracted from various wind tunnel and atmospheric 
boundary layer records, allows an accurate deduction of multiscaling exponents and cumulants. 
These exponents serve as the input for parametric estimates of the probabilistic cascade generator. 



Q 

00 

>>: 
i 

• ■ I. INTRODUCTION 

CZ2 

O 

The inertial-range dynamics of fully developed turbulent flows reveals an ubiquitous multiscale character, driven 
£>v by large-scale forcing on one end and controlled by fluid viscosity and small-scale dissipation on the other. For the 
intermediate scale range, the underlying Navier-Stokes equations do not disclose any other distinguishing length scales, 
thus suggesting a picture of a scale- invariant inertial-range dynamics This scale-invariance should reflect itself 

in the scaling of structure functions, which are moments of velocity increments constructed from the measured velocity 
data. However, measured structure functions in shear flows seem to show only an approximate multiscaling character, 
even at the largest accessible Reynolds numbers, where the large-scale and the dissipation scale are separated by 
five orders of magnitude 0, 0, 01 . The reason for this distortion of scaling seems to be the sensitivity of structure 
functions to the mean shear that is inevitably present in most natural flows at high Reynolds numbers. While schemes 
have been proposed to account for these effects [g, , it has been demonstrated recently that the lowest-order 
two-point correlation function of the surrogate energy dissipation of a high-Reynolds-number atmospheric boundary 
layer reveals a far more convincing power-law scaling over the entire inertial range, without having to resort to shear 
corrections of the sort needed for structure functions. This finding suggests that the phenomenological picture of 
""^5 ' the self-similar turbulent energy cascade may contain more truth than previously anticipated. Consequently, it is 
O \ reasonable to reexamine simple empirical energy cascade models and extend to a higher level the comparison between 
* J5 ' the models and experimental data. This is the purpose of the paper. 
£>V Random multiplicative cascade processes (RMCP) present a particularly simple geometrical picture of the energy 
. cascade 0] and are a natural archetype for multiscaling. In their binary version, for example 0, they introduce 
Oh- a hierarchy of length scales lj — L/2 3 and a random multiplicative cascade generator is used to transport the 
^ , energy flux from the integral scale L through the inertial range scales L > lj > rj down to the dissipation scale 
r\ — L/2 J . Referring to the ultrametric branching structure of binary RMCPs, TV-point correlation functions of 
arbitrary order have been calculated analytically with generating function techniques 0, 0, fpH Il3| . In terms of 
the ultrametric two-point distance, which measures the number of cascade steps necessary to reach the last common 
branching, these models reveal rigorous multiscaling. However, in this form, the two-point RMCP statistics cannot be 
compared directly to quantities that can be extracted from the data because, from an experimentalist's perspective, 
the two-point correlations are expressed as functions of the Euclidean two-point distance. In this respect, RMCPs are 
incomplete and have to be supplemented by model-dependent extensions, stating the conversion of the ultrametric to 
Euclidean statistics 0,0,0. It may be expected that this unavoidable conversion will lead to some deviations from 
rigorous multiscaling of correlation functions for two-point distances within the inertial range. The challenge is to find 
the degree to which the theoretical, extended-RMCP, two-point correlation functions can match their experimental 
counterparts. 

If this comparison is satisfactory, the two-point correlation statistics of the energy dissipation will allow for an 
unambiguous extraction of multiscaling exponents, these being related to moments of the RMCP generator. In 
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addition, two-point cumulants of logarithmic energy dissipation [lEj allow us to access properties of the cascade 
generator in another way. Taken together, reliable parametric estimations of the underlying cascade generator can 
then be given. This would be the best we can do to achieve a data-driven derivation of the turbulent energy cascade 
generator and to settle the question of whether the latter is bimodal [Til ITsf , log- normal [T^. l20j , log-stable 0, |2^| , 
log-infinite divisible |23|, |2J, |25| or another type of distribution. None of those prototype cascade generators has been 
directly derived from data. 

Note that the breakup coefficients, sometimes called multipliers, have once been thought to represent a direct 
approach to derive the RMCP cascade generator from data (26L 1271 l28t I29L l30| . In fact, generator- like scale- invariant 
distributions of breakup coefficients have been observed, but conditional distributions have been found to exhibit 
scale correlations. In a series of papers |0, these findings have been fully explained within ultrametric 

hierarchical RMCPs, once those are analyzed from an experimentalist's perspective, including unavoidable small- 
scale resummation and restoration of spatial homogeneity via the ultrametric-Euclidean conversion. This work has 
demonstrated that cascade generators and distributions of breakup coefficients are not directly related. 

The structure of this paper is as follows. Section II provides basic information about the two-point correlation of the 
lowest order for the surrogate energy dissipation in four fully developed turbulent flows — one atmospheric boundary 
layer and three wind-tunnel flows. This information serves to provide the needed guidance for the subsequent RMCP 
modelling of Sec. Ill, where analytic expressions for two-point correlators and cumulants are derived, taking into 
account the ultrametric-Euclidean conversion. More information on data analysis is provided in Sec. IV, leading in 
Sec. V to further comparisons between prototype distributions and experiments. As will also be demonstrated in 
this section, the extracted properties of the generator allow us to make parametric estimations of the RMCP cascade 
generator. A conclusion and outlook is given in Sec. VI. 

II. DATA ANALYSIS I: BASIC FACTS ON TWO-POINT STATISTICS 

We analyze four different data sets, three of which have been recorded in a wind tunnel |33| and the fourth in an 
atmospheric boundary layer [34| about 35 m above the ground. We will refer to them as wl, w2, w3 and a, respectively. 
Characteristic quantities of all data sets are summarized in Tabled The Reynolds number R\ = y/ (u 2 )\/i> is based on 
the Taylor microscale A = y (u 2 ) / ((du/dx) 2 ); v is the kinematic viscosity and u is the strcamwise velocity component. 
Upon the application of the frozen flow hypothesis, the recorded time series were converted into one-dimensional spatial 
series. The energy spectra of all four records reveal more or less the typical 5/3 slope in the inertial range. In contrast 
to wind-tunnel records, that from the atmospheric boundary layer reveals a white-noise behavior at very small scales; 
this noise, which comes from detailed electronic circuitry, has been removed by an appropriate Wiener filter. The 
energy dissipation was then calculated as the surrogate amplitude e — 15^(^) 2 . Various tests were made to ensure 
that the effect of Wiener filtering were not consequential. 

Figure 1 illustrates the lowest-order two-point correlator ri,i(d — \x2~x\\) — (e(xi)e(x2))/((e(xi))(e(x2))) sampled 
from the four different experimental records. Well inside the inertial range r\ <C d <C L the two-point correlators reveal 
a power-law behavior 7*1,1 ~ (L/d) Tl ' 1 . Power-law fits are indicated by the shifted broken straight lines; see also the 
insets, where the local slopes are shown. The resulting scaling exponents are 7*1,1 = 0.15 (wl), 0.14 (w2), 0.18 (w3) 
and 0.20 (a). Note that there is a Reynolds number dependence of this exponent. This has been explored in greater 
detail elsewhere [35|. 

For the records with the largest Reynolds number, there is a large scale range for which the two-point correlator 
exhibits a rigorous power-law scaling. However, with decreasing Reynolds number this scale range becomes smaller. 
As a rule of thumb, we observe that a good scaling range is confined between d w 20—307? and « 0.5L. If we 
could understand precisely the deviations from the power-law scaling beyond this intermediate inertial range, a more 
satisfactory extraction of multiscaling exponents would be possible, especially for turbulent flows with moderate 
Reynolds numbers. In the small-distance d < 20— 3O17 region, however, the correlations observed exceed the power- 
law extrapolation; consult Fig. 1 again. As has been explained in Ref. this enhancement is a consequence of the 
unavoidable surrogacy of the experimentally measured energy dissipation. Without knowing the correct small-distance 
behavior based on the proper energy dissipation, a theoretical modelling of two-point correlations in the dissipative 
regime is not very meaningful. We are thus left to inspect only the remainder of the inertial range. Fortunately, this 
range is sufficiently large, especially for the atmospheric data. 

At d « L, the two-point correlator has not yet converged to unity, which is its asymptotic value as d— >oo. For all four 
data records inspected, the decorrelation length appears to be around Z/dec ~ 4L and matches the length observed 
in the autocorrelation function of the streamwise velocity component. These findings suggest that the two-point 
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correlator can be described as 

ri,i(20-30»j<d<Ldec)=Oi,i( : ^ S ) hAd/L dcc ) . (1) 

For two-point distances d much smaller than the decorrelation length the finite-size function converges to unity; that 
is, <C idoc) — ► 1- One goal of this paper is to qualitatively and quantitatively reproduce this functional form 

with an extended modelling based on random multiplicative cascade processes. This approach, which is the subject 
of the next section, naturally suggests a physical interpretation of the decorrelation length Ldec- 



III. TWO-POINT STATISTICS OF RANDOM MULTIPLICATIVE CASCADE PROCESSES 
A. Binary random multiplicative cascade process 

In its simplest form an RMCP employs a binary hierarchy of length scales lj = L casc /2 J . In the first cascade step 
the parent interval of length Iq = L caac is split into left and right daughter intervals, both of length ly. In subsequent 
cascade steps, each interval of generation < j < J — 1 is again split into a left and right subinterval of length 
lj + i = Once the dissipation scale r] — L casc /2 J is reached, the interval splitting stops and has resulted into 2 J 

spatially ordered intervals of smallest size r\. It is convenient to label them as well as their ancestors according to the 
binary notation k^' = ki, k%, . . . , Kj. The label refers to the hierarchical position of an interval of generation j, where 
Ki = or 1 stands for the left or right interval, respectively. 

The binary interval splittings go together with a probabilistic evolution of the energy-flux field. From generation j 
to j + 1 the field amplitudes propagate locally as 

n re «),i = 9«ti),in K o) . (2) 

The two random multiplicative weights gi ft = Q K U) o and bright = Q K u),ii with mean (qieh) = (bright) = 1, are drawn 
from a scale-independent bivariate probability density function p(qi e { t , bright), which is called the cascade generator. 
Initially, corresponding to j = 0, the iteration (J2J) starts with a given large-scale energy flux II, which might itself be a 
random variable fluctuating around its normalized mean (II) = 1. After the last iteration J — 1 — > J, the energy-flux 
amplitude is interpreted as the amplitude £ K {j) — II K (./) of the energy dissipation supported at the interval of length 
rj at position n^. As a result of ©, this amplitude is a multiplicative sum of the random weights, given by 

£ki,...,kj Qki,...,kjQki,...,kj-i ' ' ' ^rei-fl ■ (3) 

In the following we assume the cascade generator to be of the factorized form p(q\ c f t , (fright) = p('Zicft)p(9right) with 
identical statistics for the left/right variable. Of course, the factorization is not the most general ansatz but, as 
already pointed out in Ref. 32], it represents a reasonable approximation: the turbulent energy cascade takes place 
in three spatial dimensions and calls for a three-dimensional RMCP modelling, respecting energy conservation. Since 
the measured temporal records come from one-dimensional cuts, the three-dimensional RMCP has to be observed in 
unity subdimension. Because of this, the RMCP appears to be non-conservative and the two multiplicative weights 
appear to be almost decorrelated and independent of each other. 



B. Two-point correlators 

Expressions for iV-point moments (e(xi) ni ■ ■ ■ e(xjsr) nN ) are easily found. They can be calculated either by a 
straightforward approach or, more formally, by an iterative construction of the respective multivariate characteristic 
function 0, 0] ; a third and more elegant approach [3(| makes use of the full analytic solution of the multivariate 
characteristic function for logarithmic cascade-field amplitudes 0, ^| . We simply state the results up to two-point 
correlations: 

(^D = <<r) J <n ni ) , (4) 

(e(x 1 ) ni e{x 2 ) n2 ) = ( q n ^) J - D ( q '^) D ( q '^) D (n ni+n2 ) . (5) 

Here, the binary notation k^ j ^ has been transformed into a spatial bin label x = 1 + Ylj—i K j^ J ^ which runs over 
1 < x < 2 J in units of r/. Two bins x\ = (k%, . . . , k.j-d, kj-d+i, ■ ■ ■ , Kj) and x 2 = (k\, . . . , kj-d, k 'j~d+ii ■ ■ ■ > K 'j) 
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are assigned an ultrametric distance D once the first J — D k's are identical and kj-d+i 7^ k 'j-d+v ^ n other words, 
after J — D common branches along the binary tree, the two bins separate into different branches. 
For the extraction of scaling exponents 

r n =log 2 (g n ), (6) 

it is enough to consider the two-point statistics (JSJ). In normalized form, the two-point correlators are found to scale 
perfectly as 



(e(x 1 ) n ^s(x 2 )^) 

(n™ i+ ™ 2 ) / (q ni+n3 ) 



(n»i)(n"2) V<<7" 1 ) (<f' 
/n«i+«2) (T 



J-D 



(n«i)(n«2) \ i D J ' ^ 

where Id = 2 D ~ lr q represents the characteristic two-bin distance corresponding to the ultrametric distance D > and 

Tn 1 ,n 2 = Tni+na ~ T rii — T n 2 ■ (8) 

For an experimentalist, the expression J7J) does not present an observable result. Different pairs of bins, all having 
an identical Euclidean distance rj<d<L casc , do not have an unequivocal ultrametric distance. Depending on their 
position within the binary ultrametric cascade tree, the two bins might share a cascade history that is long (small D) 
or short (large D). Consequently, as an experimentalist analyzes the two-point statistics in terms of d, the ultrametric 
expression Q has to be averaged over all D that contribute to the same value of d. In order to perform this conversion 
from an ultrametric to an Euclidean distance and, by this means, to restore spatial homogeneity, we introduce the 
discrete conditional probability distribution 






(KD<A=\log 2 d\) 


1 - d2- D 


(D=A) 


d2- D 


(A<D<J) 





(J<L><oo) 


62~ J 


(L>=oo) 



p(D\d) = { d2-° (A<D<J) (9) 



of finding the ultrametric distance D for a given Euclidean distance d in units of r\ |15| . This expression has been 
derived by employing the chain picture of independent cascade configurations; consult Fig. 2. It roughly goes as 
p(D\d) ~ 2 log2 d ~ D . The sum ^2 D=0 p(D\d) = 1 — d2~ J does not add up to unity, since p(oo\d) = d2~ J represents the 
probability that the two bins belong to different L casc -domains. 

Since the one-point statistics (e(x) n ) = (e(x + d) n ) do not depend on the spatial index x, the ultrametric-Euclidean 
conversion of the normalized two-point correlator Q leads to 



a (d^0) = Vp(£>|d)r ni , na (D)+p(oo|d) 




(10) 

This expression holds for every rj<d<L casc . For d=0 and d>L casc the normalized two-point density simply becomes 
r nun2 (d=0) = (n n i +n >)(q n i +n *) J /((n"i)(<Z ni ) J (II n *)(0 J ) and r ni . n2 {d>L casc ) = 1, respectively. The two-point 
density (|1U|) does not reveal perfect scaling anymore. Usually the second term, scaling as d/L casc , is small when 
compared to the first term except for d~L casc . The modulations, observed for the first term, are an artifact of the 
discrete scale invariance |37j of the binary random multiplicative cascade model implementation. In the following we 
will discard these modulations by first considering only dyadic distances d m = L casc /2 m with integer 0<m< J, and 
then switching again to continuous d by interpolating between the discrete d m . The expression (|1(JH then simplifies to 

/icasc\ T " 1 '" 2 

Tni, 712(d) — ^ni,n2 ( ~\ I /ni,ri2 (d/-^casc) ; (H) 
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with 



and the finite-size scaling function 



/TT"i+™2\ / I n ni+n 2 \ 

I 1 2 } q (12) 



d 



1+Tm, 



f ni ,nML C3S o) = 1 + «,„ 2 - 1) ( — J (13) 



Figure 3 compares the expressions ()l()|l and (jl 1|> for the order rii = n<x = 1. 

The finite-size scaling function has the property f ni ,n 2 (d<^L casc ) = 1 as long as the condition 1 + r ni .„ 2 > or, 
equivalently, (q ni+n2 )/((q ni )(q n2 )) > 1/2 is fulfilled. This is the case for all positive combinations ri\ > 0,n 2 > 0. 
However, combinations with negative orders do exist, for which the second term on the right hand side of l|13() 
then dominates over the first term in the limit d/L casc — > 0. This implies that the normalized two-point density 111(1 
asymptotically scales as r niv1l2 (d) ~ (L casc /d) _1 , giving rise to the effective scaling exponents n = sup{ — 1, T ni .„ 2 }. 
This scaling transition is again a pure consequence of the ultrametric-Euclidean conversion. More discussions on this 
subject can be found in Refs. pL4L I38L l39| | . 

Upon studying the expression (|13fl more closely, we realize that two effects, ultrametric-Euclidean conversion and 
large-scale fluctuations, contribute to the finite-size scaling function. They have a tendency to cancel each other. 
Once we have 

/TT™i+ n 2\ / n ni+n 2 \ 

n Vl ' - 1 , (14) 



(n n i)(n™2) (q ni )(q ni - 

the finite-size scaling function becomes f ni ,n 2 (d/ L c &sc) — 1 exactly, showing no ci-dependence. 

We also wish to point out an interesting mathematical observation following from the specific expressions (|ll|) - (|13fl . 
Since the finite-size scaling function l|13|) reveals the simple scaling behavior 

(f ni ,n 2 (d/L casc ) - 1) (^V" 1 '" 2 - (j^-X (15) 



with [i — 1, we find 



r n u n 2 (d) ~ -^r nun2 (£_d) = <X nun2 ^1 - — j ( j - I Ki) 



where the normalized two-point correlator with the rescaled two-point distance £d has been subtracted from itself. 
As a function of the two-point distance d this quantity exhibits a rigorous power-law behavior with scaling exponents 
T nii „ 2 , which is independent of the chosen rescaling parameter £ and is free of large-scale effects. 



C. Two-point cumulants 

Because experimental data yield limited statistics, the two-point correlation densities Hll|) will be restricted to 
lowest orders 1 < ri\ +ri2 < 3 or 4. This limits us to indirect information on the cascade generator, namely the scaling 
exponents t\, T2, T3 and, perhaps, T4 of JHJ ■ In order to do better, we need to accumulate additional and complementary 
information. In fact, as proposed already in Ref. |15| . this can be achieved by switching to the logarithmic amplitude 
e(x) — > lne(a;), and from two-point correlation densities to two-point cumulants 

Cis(x 2 -x 1 ) = (\n£(x 1 )\ne(x 2 )) - (\ne(x)) 2 , 



C2,i(x 2 — xi) = (in 2 e(xi) lne(x 2 )) — 2 (lne(a;i ) lne(x 2 )) (\ne(x)} - (in 2 e(x)) (\ne(x)} + 2 (lne(x)) J 

Qn 1 +n 2 



Cn una {x2-xi) = Q Xnidx , l2 ^ (e{ Xl ) Xl e(x 2 ) X2 ) 



(17) 



Ai=A 2 =0 



Explicit RMCP expressions have already been derived in Ref. [l5( within the ultrametric view as well as the converted 
ultrametric-Euclidean view. Here, we summarize only the latter lowest-order results, which hold for 7/ < d < £ casc : 

C n ^ hl (d)=G 1 (J,d)c n + G (J,d)(\n n U) c . (18) 
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The geometric functions G n (J, d) = (1/2 J ) X)d=i(^ — D) n p(D\d) are related to moments of the conditional probability 
distribution © and arc fingerprints of the hierarchical RMCP tree structure. They are given by the expressions 

G (J,d) = 1- 



L 



case 



G 1 (J,d) = (J -A) -2 d - (2- A -2- J ) « log 2 fc)-2 + 2-^-, (19) 



J CilSC 



with the last step neglecting small log-oscillations. The cumulants of the logarithmic multiplicative weight 

d n Q(X) 



z n = (ln n g) c = 



(20) 

A=0 



arc generated by the logarithm of the Mellin transform of the cascade generator, i.e., 

Q(A)=m(| dqp(q)q X ^j . (21) 
The cumulants (In™ II) c of the initial large-scale fluctuation are defined analogous to (In™ q) c . 

D. Multifractal sum rules 

The cumulants c„ of 1)20(1 and the scaling exponents r„ of © are not independent of each other. Combining Eqs. 
©, and <ED, we arrive at 

Q(n) =ln2r„ = ^ Cfc ^T ' ( 22 ) 
fc=i 



In the lowest order, this translates to 



ln2ri=0 = Cl + ^ + 5 + ... . (23) 
2 b 



These multifractal sum rules can be used, for example, to estimate ci, which can not be extracted from the two-point 
statistics <|18ll . 

Another approach to estimate the value of c\ is given by the well-known replica trick: 



n , dQ(X) 
ci = (Inq) - 



d(q X ) 



\=0 dX 



Another form of sum rules follows from (|22|l and states that 



= Um {qX) 1 = lim — . (24) 

A^O A A^O A 

A=0 



It can be seen as a generalization of (|24|l . 



n k n 2 



E n ii 
° k+1 ~kl =Cl+ nC ' 2 + ~Y C3 + ■■■ ■ ( 5 ' 



IV. DATA ANALYSIS II: MORE ON TWO-POINT STATISTICS 

In this section, we provide further aspects of data analysis that are relevant. The goal is three-fold: first and 
already pointed out in Sect. II, to test the RMCP expression 1(11(1 with the proposed finite-size scaling for two-point 
correlators, second, to test the expression ((18(1 for two-point cumulants derived from the RMCP theory, and, third, 
to extract reliable values for the scaling exponents r„ and cumulants c n from various turbulent records discussed in 
Sec. II. 

The expressions ((llll and ((18(1 come with parameters L casc , r niiTl2 , a„ 1>Jl2 , c„ and (ln"II) c . The parameter £ casc 
depends on neither the order n\, n 2 nor the choice of the two-point statistics, i.e., whether we use the correlator or the 
cumulant. The fits of ((11|) and ((18(1 to their counterparts from experimental records have to respect this independence. 
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In addition to the common parameter L casc , the fit of each order has two more parameters: these are T ni;n2 , o, ni n2 for 
the two-point correlator and c„, (In™ II) c for the two-point cumulant. We demonstrate the quality of fits by choosing 
the data record w2. For this data set, two-point correlators with orders from nx,ri2 = 1,1 to 2, 2 are illustrated in 
Fig. 4, while the two-point cumulants with orders ni,n 2 — 1,1 and 2,1 are illustrated in Fig. 5. Except for very 
small two-point distances r\ < d < 20— 30?7, where, as already noted, the surrogacy effect of the energy dissipation 
corrupts the experimental two-point statistics 0] , the agreement between the experimental two-point correlators and 
cumulants and the best-fit expressions and i|18fl is remarkable. Table II lists the best- fit parameters L casc , r„ and 
c n . Note that, due to © and n = 0, the scaling exponents T nijn2 have been converted into r 2 = T\ t i, T3 = ti,\ + t 2 
and r 4 = T34 + r 3 = t 2j2 + 2r 2 . For w2 the two values r 4 = 0.79 and 0.77, the first value extracted from 7*3, 
and the second from r 2i2 (<i), are consistent with each other, although the statistical convergence of the two two-point 
correlators of order n\ + n 2 = 4 is already beyond the limit of acceptability. For the record w3 a similar statement 
can be made, but the other records wl and a are definitely confined to m + n 2 < 3. Their best-fit parameters £ C asc, 
T n and c n are also listed in Table II. 

The cumulant c\ cannot be directly extracted from two-point cumulants. However, an indirect extraction is possible 
via the relations Ij23(l and (|24|l . Using the numerical values of c 2 and C3 determined already, the truncated multifractal 
sum rule (12; jf) leads to the C\ values listed in the second-last column of Table 2. Input for the replica-trick formula 124|) 
are the scaling exponents T\, r 2 and T3 determined already. The application of a cubic spline for the extrapolation 
results in the last column of Table 2. If a linear spline is used instead, the resulting c\ values of the records wl, w2, w3 
shift closer to the values of the second-last column, but would do worse for the record a. For an order-of-magnitudc 
estimate, however, it is safe to say that c\ ~ 0.05 for the records wl, w2, w3 and c\ w 0.07 for the record a. 

This comparison of the RMCP theory with data demonstrates that the parameter L casc is a meaningful quantity and 
deserves some fundamental consideration. For two-point distances 20 — 3O77 < d < L casc , the fitted expressions and 
111 81 are in good qualitative agreement with their experimental counterparts. At d — L casc all two-point correlations 
decorrelate and become identical to their asymptotic values corresponding to d — > 00; L casc > L is somewhat larger 
than the operationally defined integral length L of Table I, which is calculated as the correlation length of the velocity 
autocorrelation function. The physical interpretation of the extracted parameter L caac is that of a turbulent cascade 
length which, according to Fig. 2, describes the spatial extension of a hierarchical RMCP domain. 

The finite-size scaling f ni ,n 2 (d/L casc ) of two-point correlators, predicted by RMCP, now allows for an unambiguous 
derivation of scaling exponents, even for fully developed turbulent flows with a moderate Reynolds number. Conse- 
quently, as a closer inspection of Table II shows, reliable statements that the intermittency exponent t 2 might show 
a weak dependence on the Reynolds number appear to be within reach. Of course, to make this statement solid, the 
analysis of many more records is needed; this effort will be reported in |35| . 



If the scaling exponents r„, or the cumulants c n , exist and are known for all orders 1 < n < 00, the binary RMCP 
generator could in principle be reconstructed via the inverse transform of (|21fl • Unfortunately, as we have seen in the 
previous section, reliable information is limited to the lowest orders. Hence, the best we can do is to use sophisticated 
parametric estimates. Section V.A lists some of the most popular parametrizations and compares their performance 
with the results listed in Table II. Section V.B introduces the so-called log-normal inverse Gaussian distribution 
|40l l4l| . which represents a broader and more flexible parametrization class, with the purpose of finding a suitable 
approximation to the true cascade generator. 



Here we list a number of popular generators p(q) for binary random multiplicative cascade processes. They all have 
the property that the expectation value (q) = 1. 
The log-normal distribution 



is classic |l9j,|20J. Its log-stable generalization has also been considered |2lU22|. but does not qualify for our purposes, 
since the cumulants c„ do not exist for this distribution beyond some order. For comparison, we will also employ the 



V. 



PARAMETRIC ESTIMATION OF THE RMCP GENERATOR 



A. Dictionary of prototype cascade generators 




(26) 
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rescaled gamma distribution 



and the asymmetric beta distribution [3/ 



Pgammate) = f^j^ ^ (27) 



The bimodal distribution 



PUph.(«) = — X—« (9- (l-ai)) + —x— «(«-(! + <»)) ' < 29 ) 
OL\ + a 2 ai + a 2 

although discrete, has also been used extensively 0, ■ Another popular discrete representative is the log-Poisson 
distribution 

(g) = E 2 ~^ ln2 )% ( g - 2^"^f) , (30) 



PlogPoisson 



which was originally derived with v\ — 2 and = 2/3 from some plausible reasoning on the structure of the most 
singular objects in fully developed turbulent flows (23L I24L l25j . 

For all parametrizations l|26() - (|30|l it is straightforward to determine analytic expressions for the scaling exponents 
and cumulants via l|2(J[l - H22[l . The free parameter of the one-parametric distributions (|26[l - i|28[) is then fixed to 
reproduce the observed intermittency exponent r 2 , listed in Table II. The two-parametric distributions (|29|l and l|3U|l 
need also conform to T3 in addition to r 2 . No further freedom is left for the scaling exponents of higher order and 
cumulants of all orders. Table III summarizes their predicted values. 

It is difficult to rate the prototype cascade generators because of ambiguity inherent in the data. Within the one- 
parametric distributions l|2t)[l - H28|) the log-normal distribution performs better: for all the records, the predicted values 
for ci and c 2 are close to the observed cumulants. However, the log-normal distribution without skewness is unable 
to reproduce the observed positive values for C3 and for the record w2 it also overestimates the scaling exponent T4. 
Furthermore, the difficulties of the log-normal distribution for high-order moments is now well known 2|. Compared to 
the log-normal distribution, the rescaled gamma distribution and the asymmetric beta distribution have the tendency 
to overestimate the first two cumulants. Furthermore, C3 is predicted with an opposite sign. Rather surprisingly, the 
simple two-parametric bimodal distribution 129(1 shows the closest agreement for all records. The scaling exponent 
T4, if observed, as well as the cumulants c\ and c 2 almost match their observed counterparts. Moreover, C3 comes 
with the correct sign, although it is about a factor 2 too low for the records wl, w2, w3 and roughly a factor 2 too 
large for the atmospheric boundary layer record. Like the distributions (|2"?|) and the two-parametric log-Poisson 
distribution overestimates the second cumulant c 2 and, except for record w3, predicts C3 with the wrong sign. It is 
interesting to note that the parameter-free log-Poisson distribution j2j with v\ = 2 and j/ 2 = 2/3 matches well the 
scaling exponents T3 and T4 of record a with the largest Reynolds number, but disagrees with all cumulants. 



B. Log- normal inverse Gaussian distribution 

A broader and more flexible parametrization class is the so-called normal inverse Gaussian distribution |4fJ. l4l| 

p(x;a,/3,»,5) = a(a,p, f M,6)s( 3 ^) K x (sas( e? x , (31) 



with s(x) = \J {\ + x 2 ), a(a, /3,/x, 5) — n~ 1 aexp(5y / a 2 — f3 2 — /3/x) and — oo<a;<oo. K\ is the modified Bessel function 
of the third kind and index 1. The domain of variation of the four parameters is given by /i G R, 6 G R + and < 
|/9| < a. The distribution is denoted by NIG(a, (3, /i, S), and its cumulant generating function Q(A; a, (3, /i, 5) = ln(e Ax ) 
has the simple form 

Q(A; a, 0, fi,S)=S [^a 2 - (3 2 - ^a 2 - (/? + A) 2 ) + /iA . (32) 

If X\, . . . ,x m are independent normal inverse Gaussian random variables with common parameters a and (3 but 
individual location-scale parameters fii and Si (i = 1, . . . , m), then x+ = x\ + . . . + x m is again distributed according 
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to a normal inverse Gaussian law with parameters a, 0, /Lt_|_ and Furthermore, we note that the NIG distribution 
(|31[l has semi-stretched tails 

p(x; a, (3, (m, S) ~ 3 ^ 2 exp (— a \x\ + [3x) (33) 

as x — > ±oo. This result follows from the asymptotic relation K v (x— >oo) ~ ^tt /2x~ 1 / 2 e~ x . 
For our purposes we assume the random multiplicative weight to be distributed according to 

lng~NIG(a,/3,£,M) , (34) 

which turns normal inverse Gaussian statistics into log-normal inverse Gaussian statistics. With (|2U[) . (1221) and H32[) . 
the scaling exponents and cumulants yield 



r n In 2 = Q(n; a, (3, fi, S) (35) 



and 



y'l — p 2 ayX — p) 1 cr(l — p z ) t >' z 

where p = (3/ a. 

For each of the records wl, w2, w3 and a, the four NIG parameters a, f3, S and p are determined so as to reproduce 
T\ = and the observed values for t-i , T3 and C2 listed in Table [H] Since the respective expressions (|35f) and 1301) 
are nonlinear, real solutions for the parameters are not guaranteed. Where complex-valued solutions resulted in the 
first attempt, which occurred for w2, w3 and a, the values for C2, T3 and T2 are relaxed, in this order and to some 
small extent, until real- valued parameter solutions are obtained. The outcome is listed in Table HVI The results are 
very close to the log-normal values listed in IIIII The log-normal inverse Gaussian distribution has the tendency to 
overestimate the fourth-order scaling exponent T4. The magnitude of the third cumulant C3 is strongly underestimated, 
so that its predicted sign shows only random scatter. Figure [H] compares the extracted log-normal inverse Gaussian 
distributions of Table II VI with the extracted log-normal distributions of Table 11111 For all four records the extracted 
distributions are very similar. 

As a summary of this section, we reiterate that the bimodal distribution (|29|l produces the best overall agreement 
with the observed scaling exponents and cumulants. However, the true cascade generator will not be discrete. From 
the set of continuous generator representatives tested, the log-normal and log-normal inverse Gaussian distributions 
perform best and about equally well. 



VI. CONCLUSION 



Random multiplicative cascade processes are able to describe the observed two-point correlation structure of the 
surrogate energy dissipation of fully developed turbulent flows beyond simple power-law scaling. Keeping in mind the 
need for a satisfactory comparison between modelling and experimental data, a useful transformation has been intro- 
duced: this transformation converts model-dependent and unobservable ultrametric two-point distances to Euclidean 
two-point distances, the latter reflecting the horizontally sampled n-point statistics of the experimental records. The 
predictions of RMCP for finite-size scaling of two-point correlation functions are confirmed by experimental data 
from three (wind tunnel) shear flows and one atmospheric boundary layer; a new physical length scale characterizing 
the upper end of the inertial range, called here the turbulent cascade length, has been shown to appear naturally. 
Furthermore, the quantitative classification of the deviations from a rigorous scaling of two-point correlators allows 
for an unambiguous extraction of multiscaling exponents, even for flows with moderate Reynolds numbers. When 
complemented with additional information extracted from two-point cumulants of the logarithmic energy dissipation, 
reliable parametric estimates of the RMCP generator have become feasible. 

RMCPs produce a consistent geometrical modelling of the self-similar turbulent energy cascade. This is further 
supported by recent investigations on the scaling part of three-point statistics [4^ | and previous investigations on scale 
correlations [HI I31H32L l43| . Thus, the self-similar and RMCP-like cascade process appears to be universal. However, 
as shown by variations among the four records considered, the strength of the cascade generator appears to depend 
on the Reynolds number and perhaps also the flow geometry. This needs further clarification by means of a separate 
and extended effort [35|; needless to say, there is a strong need for the experimentalists to produce clean and longer 
records with converged statistics. 
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data set 


Rx 


L/rj 


-^record / L 


\/ V 


Ax/rj 


wl 


306 


484 


102500 


35 


1.97 


w2 


493 


968 


193500 


44 


2.79 


w3 


1045 


2564 


77500 


64 


2.97 


a 


9000 


5xl0 4 


1000 


187 


1.29 



TABLE I: Taylor-scale based Reynolds number R\, integral length scale L in units of the dissipation scale r?, record length 
irccord, Taylor microscale A and resolution scale Aa; of three wind tunnel (wl, w2, w3) |33l and one atmospheric boundary 
layer (a) |34| records. 



data set 


^casc J f] 


T2 




T4 


C2 


C3 


ciJ23 


dim) 


wl 


1873 


0.15 


0.46 




0.100 


0.042 


-0.057 


-0.045 


w2 


3069 


0.15 


0.42 


0.78 


0.095 


0.045 


-0.055 


-0.046 


w3 


7117 


0.17 


0.52 


0.98 


0.099 


0.060 


-0.059 


-0.046 


a 


322500 


0.21 


0.58 




0.149 


0.015 


-0.077 


-0.068 



TABLE II: Parameter values resulting from least-square fits with expressions (111! and 11811 . L c 
set. The last two columns represent the estimates 12311 and 12 It for the cumulant c\. 



has been fixed for each data 



10 u 




10 u 



10" 



10" 



10' 
d/ii 



10" 



1(f 



10" 




10' 



10" 



10 



d/n 



10 



10" 




10' 



d/r| 



10' 




10 



d/n 



10 



10" 



FIG. 1: Two-point correlator {e(x + d)e(x)} / (e(x)} 2 of the surrogate energy dissipation extracted from the records wl, w2, w3 
and a. Power-law fits are indicated by the shifted broken straight lines. The inset figures illustrate the local slope. 
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log-normal distribution (l-Mil 



data set 


a 


T3 


T4 


Cl 


C2 


C3 


wl 


0.33 


0.46 


0.93 


-0.054 


0.107 


0.000 


w2 


0.32 


0.44 


0.87 


-0.050 


0.101 


0.000 


w3 


0.34 


0.51 


1.03 


-0.059 


0.119 


0.000 


a 


0.38 


0.63 


1.26 


-0.073 


0.145 


0.000 


gamma distribution 1271 


data set 


7 


T3 


T4 


Cl 


C2 


C3 


wl 


8.84 


0.45 


0.87 


-0.058 


0.120 


-0.014 


w2 


9.41 


0.42 


0.82 


-0.054 


0.112 


-0.013 


w3 


7.94 


0.50 


0.96 


-0.064 


0.134 


-0.018 


a 


6.39 


0.60 


1.16 


-0.080 


0.169 


-0.029 


beta distribution 1281 


data set 


P 


T3 


T4 


Cl 


C2 


C3 


wl 


7.61 


0.44 


0.85 


-0.059 


0.124 


-0.019 


w2 


8.11 


0.42 


0.81 


-0.055 


0.116 


-0.017 


w3 


6.82 


0.49 


0.94 


-0.066 


0.139 


-0.025 


a 


5.47 


0.60 


1.13 


-0.083 


0.178 


-0.040 


bimodal distribution 1291 


data set 


ai = 


«2 


T4 


Cl 


C2 


C3 


wl 


0.22 


0.52 


0.88 


-0.049 


0.092 


0.025 


w2 


0.24 


0.44 


0.80 


-0.049 


0.094 


0.018 


w3 


0.21 


0.61 


1.00 


-0.052 


0.094 


0.033 


a 


0.31 


0.50 


1.06 


-0.075 


0.144 


0.026 


log-Poisson distribution 13011 with v\ = 


= 2, v 2 


= 2/3 




T2 


T~3 


T4 


Cl 


C2 


C3 




0.22 


0.59 


1.06 


-0.10 


0.228 


-0.092 




log-Poisson distri 


bution EOJ 




data set 


V\ 




T4 


Cl 


C2 


C3 


wl 


109.84 


0.96 


0.90 


-0.055 


0.111 


-0.004 


w2 


13.75 


0.90 


0.82 


-0.054 


0.112 


-0.012 


w3 


1548.7 


1.01 


1.03 


-0.059 


0.117 


0.001 


a 


4.61 


0.79 


1.09 


-0.085 


0.184 


-0.044 



TABLE III: Fitted parameters for a few prototype cascade generators and their predicted values for the remaining scaling 
exponents r n and cumulants c„. 



log-normal inverse Gaussian distribution 1341 



data set 


Q 


P 


S 


A 4 


T2 


T3 




Cl 


C2 


C3 


wl 


10.98 


0.94 


1.09 


-0.14 


0.15 


0.46 


0.95 


-0.051 


0.100 


0.002 


w2 


17.28 


-5.46 


1.56 


0.47 


0.15 


0.43 


0.85 


-0.052 


0.106 


-0.006 


w3 


27.26 


17.80 


1.17 


-1.06 


0.16 


0.52 


1.11 


-0.052 


0.099 


0.012 


a 


6.99 


-2.44 


0.92 


0.27 


0.20 


0.60 


1.19 


-0.076 


0.160 


-0.027 



TABLE IV: Fitted parameters of the log-normal inverse Gaussian cascade generator and their predicted values for the scaling 
exponents r n and cumulants c n . 
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10° 10 1 10 2 10 3 

d/r| 

FIG. 3: Comparison of the expressions 1101 and [11 1 1 for the order n\ = ri2 = 1, showing that the log-oscillations inherent to 
HUH remain small. Parameters have been set as follows: L C!lBC /rj = 2 10 , T2 = 0.20 and II = 1. 
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FIG. 4: Best fits of expression 11 U to two-point correlators extracted from data set w2. Correlation orders are n\,ni = 1, 1, 
2, 1, 3, 1 and 2, 2. The parameter Lease has been fitted such that it is the same for all orders of two-point correlators as well as 
two-point cumulants, with the latter illustrated in Fig. 5. The inset figures illustrate the local slope. For comparison, power-law 
fits with the extracted scaling exponents listed in Tab. II are shown as the shifted dashed straight lines. 
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1.5 




FIG. 5: Best fits of expressions 1181 to two-point cumulants extracted from data set w2. Correiation orders are ni,ri2 = 1,1 
and 2, 1. The parameter Lease has been fitted such that it is the same for ail the orders of two-point cumulants as well as 
two-point correlators, with the latter illustrated in Fig. 4. 




FIG. 6: Comparison of the extracted log-normal distribution 1261 1 (solid) and log-normal inverse Gaussian distribution 13411 
(dashed) for the records wl, w2, w3 and a. Respective parameter values have been taken from Tables fTTTI and HVI 



